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We propose a concept of a quantum cascade laser based on transitions of bosonic quasiparticles 
(excitons and exciton-polaritons) in a parabolic potential trap in a semiconductor microcavity. This 
laser would emit terahertz radiation due to bosonic stimulation of excitonic transitions. Dynamics 
of a bosonic cascade is strongly different from the dynamics of a conventional fermionic cascade 
laser. We show that populations of excitonic ladders are parity-dependent and quantized if the laser 
operates without an external terahertz cavity. 
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Quantum cascade lasers (QCLs) are based on subse- 
quent intersubband transitions of electrons or holes in a 
Wannier-Stark ladder formed in a semiconductor super- 
lattice subject to an external electric field [ll-Q- Emitted 
terahertz (THz) photons are polarized perpendicularly to 
the plane of the structure and propagate in-plane (which 
is referred to as wave-guiding or horizontal geometry). 
QCLs differ from conventional lasers as they do not re- 
quire inversion of electronic population for every partic- 
ular transition. Still, this is a fermionic laser, where the 
quasiparticles emitting light obey Fermi statistics. Re- 
cently, several proposals for bosonic THz lasers based 
on exciton-polaritons have been published 0, [B[ . These 
sources are expected to generate THz light beams propa- 
gating in the normal-to-plane direction (vertically) with- 
out external THz cavities @. The emitted radiation is 
stimulated by the final state (exciton-polariton) occupa- 
tion, which is a purely bosonic effect. 

Here we propose a concept of a bosonic cascade laser, 
which combines the advantages of QCL (emission of mul- 
tiple THz photons for each injected electron) and exciton- 
polariton lasers (no need for a THz cavity, low threshold) . 
We consider an exciton cascade formed by equidistant en- 
ergy levels of excitons confined in a parabolic trap in a 
semiconductor microcavity. Parabolic traps for exciton- 
polaritons have been experimentally demonstrated, and 
an equidistant spectrum of laterally confined exciton- 
polariton states has been observed. There are several 
ways to realize such traps including specially designed 
pillar microcaviti es Bl , strain induced Q and optically 
induced traps |9l4ll|. A particularly promising varia- 
tion of these designs would be a microcavity with a large 
parabolic quantum well embedded. We consider the weak 
coupling regime where the optical mode is resonant with 
the m th exciton level to allow efficient pumping. The 
other energy levels of the confined excitons would be un- 
coupled to the cavity mode, forming a dark cascade ideal 
for a high quantum efficiency device due to the long ra- 



diative lifetime. This device would emit radiation polar- 
ized in the direction normal to the quantum well plane 
and propagating in the cavity plane in the wave-guiding 
regime. In this Letter we formulate a kinetic theory of 
bosonic cascade lasers and calculate the matrix elements 
of terahertz transitions in realistic parabolic traps. 

The occupation numbers of exciton quantum confined 
states and the THz optical mode in our cascade laser can 
be found from the following set of kinetic equations (0 is 
the state with the lowest energy, m is the state with the 
highest energy, which is resonantly pumped): 
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Here W = Wo(tithz + 1) is the THz emission and 
W' = WoriTHz is the THz absorption rate, uthz is the 
THz mode occupation, and t^hz is the THz mode life- 
time. We assume that the matrix element of THz tran- 
sition is non-zero only for neighbouring stairs of the cas- 
cade and that it is the same for all neighbouring pairs. 
This simplifying assumption allows for the analytical so- 
lution of Eqs. (JT|) — ((3j ■ It can be easily relaxed in the 
numerical calculation, accounting for the specific condi- 
tions of particular experimental systems. 
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Wc first consider the case where there is no THz cavity, 
and assume that THz photons leave the system immedi- 
ately such that 71thz = 0. The solid curves in Fig.Q]show 
the dependence of the mode occupations on the pump 
power in this case, which were calculated by numerical 
solution of equations Jl]) - §5§ for the steady state. For 
increasing pump power we see that the modes become 
occupied one-by-one and a series of steps appears, each 
corresponding to the occupation of an additional mode. 
In the limit, Wqt -C 1, the position of the steps is given 
by P/Wq = n 2 / (Hot) 2 where n is a half-integer. For high 
pump powers, where all modes are occupied two different 
behaviours of the modes can be identified: 0th, 2nd, 4th 
etc modes continue to increase their occupation with in- 
creasing pump power, while 1st, 3rd, 5th etc modes have 
a limited occupation. This effect persists independently 
of whether an even or odd number of modes is considered 
in the system. 
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FIG. 1: (color online) Dependence of the mode occupations 
in the absence of a THz cavity on pump intensity, calculated 
numerically from the kinetic Eqs. (HJ — (EJ (solid curves) and 
analytically from Eqs. ([SJ, (0 and (J7J) (dashed curves). The 
vertical gray lines correspond to the step locations given by 
P/Wo = n 2 /(W T) 2 where n is a half-integer. Parameters: 
Wot = 3 x 1CT 6 , m = 15, n THz = 0. 



Qualitatively, our results can be understood as follows. 
Every mode in the chain experiences both a gain and a 
loss. The last mode in the chain is unique since it only 
experiences loss due to the finite lifetime rather than THz 
emission. Since it experiences loss only due to the life- 
time, we can expect that the last mode is strongly oc- 
cupied in the limit of high pump power. This means 
that the second-to-last mode experiences a strong loss 
due to stimulated scattering to a highly occupied state. 
Thus the second-to-last mode has a much smaller occu- 
pation. The third-to-last mode then experiences only a 
small loss due to stimulated scattering and so can again 
have a large occupation. The series repeats such that 
alternate modes have high and low intensity, with the 
highly occupied modes introducing a fast loss rate that 
limits the occupation of low intensity modes. 

Quantitatively, Eqs. ([J) - (j3J) can be solved analyti- 
cally in the steady state, where J\ Nk = Pt, under the 



assumption that Wqt <C 1, N m 1 and titHz = 0. In 
this regime, Nq depends linearly on the pump intensity: 
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where |~n] denotes rounding up to the nearest integer. 
Our approximation makes sense when Nq is positive, i.e., 
when Pt > m 2 / (2Wot) 2 . The populations of modes with 
even index also depend linearly on the pump intensity: 
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The populations of the odd modes are: 
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Results from Eqs. |(5Jl, © and ([7]) arc compared to the 
numerical results in Fig. [T] It is also possible to write an 
equation for the THz emission rate: 
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Results from Eq. ([9]) arc compared to numerical calcula- 
tion of the THz emission rate in Fig. @a. For high pump 
powers, the rate increases linearly with the pump power. 
This represents a limit to the quantum efficiency of THz 
emission, which is given by the THz emission rate divided 
by the pump rate, P, and plotted in Fig. ^2)p. While the 
presence of the bosonic cascade allows quantum efficien- 
cies exceeding unity, the quantum efficiency is limited 
under high pump powers to [m/2] , for Wot <C 1. 

Figure [2J3 shows the typical time dependence of the 
mode occupations after the pump is switched on (as- 
sumed instantaneously). The enhancement of the scat- 
tering via stimulated processes allows the system to reach 
equilibrium in a time less than 1/Wb- The presence of 
multiple, dynamically changing effective scattering rates 
causes an initially chaotic dynamics. 

In the presence of a THz cavity, macroscopic numbers 
of THz cavity photons allow further stimulated enhance- 
ment of scattering between the modes. In this case, where 
upward transitions arc allowed in addition to downward 
ones, the steps observed in Fig. [1] are washed out as shown 
in Fig. [3] In contrast to Fig. [TJ all modes continue to 
increase their population with increasing pump power. 
Figure [3Ji shows the THz emission rate as a function 
of the THz photon lifetime. For fixed pump intensity, 
higher THz emission rates are observed than in Fig. [2ji. 
Figure HJd shows the dependence of the THz emission 
quantum efficiency on the THz photon lifetime. Increas- 
ing of the THz photon lifetime is able to increase the 
quantum efficiency, however, not beyond the limit of m/2 
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FIG. 2: (color online) a) Dependence of the THz emission 
rate on pump intensity in the absence of a cavity, for different 
numbers of modes in the chain (values of m are marked on 
the plot). Solid curves show results from numerical solution 
of Eqs. (TTJ) - ; Dashed curves show the results of Eq. iJSJ, 
valid for high pump powers, b) Dependence of the quantum 
efficiency on pump intensity (the values of m are the same 
as in (a)), c) Time dynamics for m = 6, P = 3 x 10 12 W. 
Parameters: Wot = 3 x 10~ 6 , wthz = 0. 
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FIG. 3: (color online) Same as in Fig. [T] with a THz cavity. 

TTHz = T/1000. 



already observed in the absence of a THz cavity for high 
pump powers. 

Let us now demonstrate the feasibility of THz transi- 
tions between the size-quantized levels of excitons in a 
parabolic quantum well. The two-particle Hamiltonian 
for an electron and a hole reads 



2fi dp 2 
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+ V(z e ,Z h ). (10) 



Here M = m e + ruh, p = m e m,h/M is the reduced mass 
of the electron-hole pair, p — r e — rh = (x, y, z) is the rel- 
ative coordinate and R = (niere+nihr^/M = (X, Y, Z) 
is the center of mass wave vector, e is the background di- 
electric constant. The quantum well potential V(z e ,Zf l ) 
is written in form: 
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FIG. 4: (color online) a) Dependence of the THz emission rate 
in the presence of a THz cavity on the THz photon lifetime. 
Different curves correspond to different pump powers (marked 
on the plot), b) Dependence of the quantum efficiency on THz 
photon lifetime for the same pump powers in (a). Parameters: 
m= 15, W r = 3x nn 6 . 



where constants A e ,h denote corresponding stiffness. It 
is convenient to rewrite V(z e ,Zh) as a function of the 
center of mass Z and relative motion z coordinates with 
the result V = V (Z, z)+Vi(Z, z), where V (Z, z) = {A e + 



A h )Z 2 + (m 2 e + m 2 )z 2 /M 2 
and relative motion and 



does not mix center of mass 



Vi(Z, z) = — (m h A e 



,A h )Zz, 
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mixes internal and center of mass degrees of freedom pro- 
vided that m,hA e ^ m e Ah- 

In what follows it is assumed that the potential is 
weak enough, hence, in the zeroth approximation center 
of mass can be quantized independently and th order 
wavefunctions have the form: 



^n, m ,hi r e,r h ) = R n i m {p)Yi m (-d,(p)F h (Z), 
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where Fh(Z) are the eigenfunctions of the potential 
(A e + Ah)Z 2 , h enumerates levels in this harmonic poten- 
tial, Rnim(p) are the 3D-hydrogen-like radial functions 
and Yi m ("d, ip) are (3D) the angular harmonics of relative 
motion, n, I and m are the quantum numbers of relative 
motion of an electron and a hole. Equation (|12|) holds if 
Kl <C 7Z, where f2 is the eigenfrequency of the center of 
mass in the parabolic potential, MCI 2 /2 = A e + Ah, and 
1Z = ixe A / (2e 2 h 2 ) is the exciton Rydberg. 

It is the mixing of the exciton states caused by per- 
turbation Eq. (fTTj) that gives rise to the THz transitions 
between the neighbouring size-quantized states in a quan- 
tum well. Making use of matrix elements of coordinate 
for harmonic oscillator and for the hydrogen atom 12|, [13 1 
we obtain for the mixing matrix elements (bra- and ket- 
denote states as, e.g. |n, I, m, h)) 

(l,0,0,/i-l|Vi|n,l,0,/i) = (1,0,0, h\Vi\n, 1,0,^-1) = 
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where ao = ^f- is the exciton Bohr radius and 

_ 2 4 (n- l)»-5/2„7/2 
jn ~ (n + l)«+5/2 • 

Using first order perturbation theory, the THz transition 
dipole moment between the neighbouring states is: 

D h = e^xlxl**) = 7T#T% (14) 
^ A /l ,±(n) V3 

where the wavefunctions of Is exciton states take the 
form in the leading order in V\ : 

*i,o,o,fc(p,ii) = i?ioo(p)>oo(^^)^(^)+ 

E A )^ n W {p)Y w {$,ip)F h±l {Z). (15) 

Here A hi ±(ri) = E ls - E np =p ftfi « £i s - £J np , and £?„ s 
/'.', : , = —IZ/n 2 . Estimation according to Eq. (fT4"]) shows 
for a = 100 A, TZ = 5 meV, Sf2 = 1 meV, m e = O.lmo, 
to/, = 0Am o and A e = Ah that 

|£) h | = eVft x 8.3 x 10 1 A, (16) 

which is comparable with 2p — > Is transition strength. 
Following Ref. 0] , this gives the transition rate for ITHz 
emission in the kinetic equations Wo ~ 5500s -1 (where 
we took the refractive index as 3 corresponding to a GaAs 
based system). For excitons with lifetime ~ 500ps, this 
gives W t ps 3 x 10~ 6 . 

Since the transition rate is proportional to |-D/i| 2 , 
Eq. (|16[) introduces a proportionality between the tran- 
sition rate and the mode indices, n, for the specific case 
of a parabolic trap with electron-hole mixing. This can 
be modelled by straightforward modification of the rate 
Eqs. (|T]) ([HI) - While this can distort the steps observed 
in Fig. [1] as shown in Fig.[5j the behaviour is similar. For 
Wqt <C 1 and m even, one can obtain analytically: 

1 m(l + m- PW T 2 )T(l/2 + k)T(m/2) 

N2k ~ AW TT{l + k)T[{3 + m)/2] ' 

(17) 

where k is an integer and r(n) is the Euler gamma func- 
tion. The other modes of the system all have the same 
intensity, N 2 k+i ~ 1/(Wot) anc ^ the THz emission rate 
is given by (2Pr - m/(W T) - N + N m ) /(W t). 

Finally, we note that the waveguiding geometry allows 
for the increase of the active volume of the laser and, 
consequently, the output power. Its drawback is a larger 
size compared to the vertical cavity laser geometry. Lat- 
eral polariton traps would allow a more compact size of 
the device, while the matrix elements of terahertz tran- 
sitions are believed to be significantly larger in parabolic 
quantum wells. 

In summary, we have proposed the use of a bosonic 
cascade to implement THz lasing with high quantum ef- 
ficiency, above unity. Such a device requires a parabolic 




FIG. 5: (color online) Same as in Fig. [T]with transition rates 
proportional to the mode index. 

trapping, which can be arranged in the plane of a semi- 
conductor microcavity via a variety of techniques or using 
parabolic wide quantum wells. Electron- hole mixing was 
shown to give rise to non-zero transition matrix elements. 
A series of steps in the mode occupations of the bosonic 
cascade was predicted for increasing pump power in the 
absence of a THz cavity. The quantum efficiency of a de- 
vice can be improved by increases in the pump power or 
the use of a THz cavity, however, the quantum efficiency 
is limited to one-half of the number of modes present. 
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